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We have studied the crossover between thermally assisted and pure quantum tunneling in single 
crystals of high spin (S = 10) uniaxial single molecule magnet Mni2-acetate using micro-Hall effect 
magnetometry. Magnetic hysteresis experiments have been used to investigate the energy levels 
that determine the magnetization reversal as a function of magnetic field and temperature. These 
experiments demonstrate that the crossover occurs in a narrow (~ 0.1 K) or broad (~ 1 K) tem- 
perature interval depending on the magnitude and direction of the applied field. For low external 
fields applied parallel to the easy axis, the energy levels that dominate the tunneling shift abruptly 
with temperature. In the presence of a transverse field and/or large longitudinal field these energy 
levels change with temperature more gradually. A comparison of our experimental results with 
model calculations of this crossover suggest that there are additional mechanisms that enhance the 
tunneling rate of low lying energy levels and broaden the crossover for small transverse fields. 
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High spin single molecule magnets Mni2-acetate 
(Mn.12) and Fes have been actively studied as model 
systems for the behavior of the mesoscopic spins [jLHL2|. 
These materials can be considered as monodisperse en- 
sembles of weakly interacting particles with net spin 
5 = 10 and strong uniaxial anisotropy. They enable stud- 
ies of both classical and quantum effects through macro- 
scopic magnetic measurements. These clusters show en- 
hanced relaxation of magnetization at regular intervals 
of magnetic field, attributed to the tunneling of magne- 
tization across the anisotropy barrier (|||. 

The temperature dependence of this process suggests 
that both classical thermal activation and quantum tun- 
neling are important || . Other significant recent experi- 
mental results on single molecule magnets include obser- 
vation of non-exponential relaxation of magnetization || 
and quantum phase interference in Fes @ ■ Spectroscopic 
experiments (EPR and inelastic neutron scattering) have 
provided important information about the magnetic en- 
ergy levels of Mni2 and Fes, relevant to understanding 
their macroscopic magnetic response |8|-[T2| . 

Experimental results, such as the series of regular steps 
in magnetic hysteresis curves of these materials, have 
been interpreted within an effective spin Hamiltonian for 
an individual cluster: 
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where the uniaxial anisotropy parameters D and B 
have been determined by EPR |l0[| and inelastic neu- 
tron spectroscopy experiments |ll[] (D = 0.548(3) K, 
B = 1.17(2) x 10~ 3 K, and g z is estimated to be 1.94(1)). 
Here TC includes terms which do not commute with 



FIG. 1. Schematic of the double well potential in Mni2. 
Solid arrows show magnetization reversal by thermally as- 
sisted tunneling, with thermal activation to an excited level 
in a metastable well and tunneling between resonant levels. 
The dashed arrow illustrates pure quantum tunneling. 



S z and produce tunneling. These mechanisms of level- 
mixing may be due to a transverse field (such as hyper- 
fine, dipolar fields or an external field, contributing terms 
like H X S X ) or higher order transverse anisotropies, for ex- 
ample C{S%+St), C = 2.2(4) x 10" 5 K Q, which is the 
lowest order term allowed by the tetragonal symmetry of 
the Mn 12 crystal. The enhanced relaxation of the magne- 
tization at certain field values is ascribed to thermally as- 
sisted tunneling (TAT) or pure quantum tunneling (QT) 
(see Fig. |TJ> . In the TAT regime the relaxation occurs 
by tunneling from thermally excited magnetic sublevels 
(rricsc = 9, 8, — 8, — 9), when two levels on the oppo- 
site sides of the anisotropy barrier are brought close to 
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resonance by the magnetic field. From the unperturbed 
Hamiltonian (1) the longitudinal (z-axis) field at which 
the levels m osc and m! become degenerate is: 

H(n, m osc ) = nH a {\ + B/D[m 2 csc + (m csc - n) 2 }}, (2) 

where n = m csc + m is the step index describing the 
bias field, and Hq = D/g z fiB is a constant (0.42 T). Di- 
rect numerical diagonalization of the spin Hamiltonian 
(1) shows that the small transverse anisotropy term does 
not significantly change these resonance fields. Note that 
a larger magnetic field is necessary to bring lower lying 
sublevels into resonance. As the temperature decreases, 
the thermal population of the excited levels is reduced, 
and these states contribute less and less to the tunnel- 
ing. Consequently, the steps in hysteresis curves shift to 
higher bias field values, and steps with larger n become 
observable. At low temperature, the tunneling from the 
ground state m GSC = 10 dominates, and the position and 
amplitude of the steps become independent of tempera- 
ture, denoted the QT regime. 

At any temperature the magnetic relaxation rate is 
dominated by the contribution of only few (one or two) 
magnetic sublevels. This is because the relaxation rate is 
the product of the thermal occupation probability and 
the tunneling rate. The probability of thermal occu- 
pation decreases exponentially with energy, while the 
tunneling probability increases exponentially approach- 
ing one at the top of the anisotropy barrier. Therefore 
the relaxation rate has a sharp maximum at the domi- 
nant level or levels. 

Recent theoretical models suggest that in the large spin 
limit a crossover between thermal activation and quan- 
tum tunneling can either occur abruptly, in a narrow 
temperature interval, or gradually, in a broader tempera- 
ture interval |l5|-|l5[] . In the first case the energy at which 
the system crosses the anisotropy barrier shifts abruptly 
with temperature from a value close to the top of the bar- 
rier to the lowest lying level in the metastable well (de- 
noted a first-order crossover). In the second case this en- 
ergy changes smoothly with temperature (second-order) 
|jl6|| . The "phase diagram" for this crossover depends 
on the form of the spin Hamiltonian, particularly the 
terms important for tunneling. In finite spin systems the 
crossover is always smeared, nevertheless the two types 
of the crossovers can be distinguished experimentally. In 
the first-order crossover there are competing maxima in 
the relaxation rate versus energy and the global maxi- 
mum shifts abruptly from one energy to the other as a 
function of temperature. In the second-order case a sin- 
gle maximum in the relaxation rate shifts continuously 
with temperature. Recent experiments have shown that 
in Mni2 the crossover occurs in a narrow temperature in- 
terval when the applied field is parallel to the easy axis of 
the sample ]l8| ]. In contrast, experiments on Fes suggest 
a second-order crossover [fill . 



In this paper we show that in Mni2 both types of the 
crossover can be experimentally observed when the ap- 
plied field has a transverse component. A transverse 
magnetic field makes the crossover more gradual and 
leads to a continuous shift in the dominant energy lev- 
els (i.e., a second-order crossover) Q. We show that 
these levels can be identified by hysteresis experiments 
and their behavior as a function of temperature can be 
studied in order to test the tunneling model described 
above. 

Our experiments have been conducted using a micro- 
Hall-effect magnetometer |2lJ in a high field 3 He system. 
Single crystals of Mni2 in the shape of parallelepipeds 
50 x 50 x 200 /im 3 were synthesized as described in 
Ref. A crystal was encapsulated in thermally con- 

ducting grease and the temperature was measured with 
a calibrated carbon thermometer a few millimeters from 
the sample. The angle 9 between the easy axis of the 
crystal and the applied magnetic field was varied by ro- 
tating the sample in a superconducting solenoid. 

In our hysteresis experiments the sample was initially 
saturated (Mo = — M s ), then the field was ramped at a 
constant rate (0.2 T/min) towards positive saturation. 
Hysteresis curves show steps and plateaus, separated by 
a field interval of approximately 0.44 T, in agreement 
with previously published results. Figure || shows a plot 
of the derivative of the magnetization dM/dH versus the 
applied field for 9 = 10° measured at different tempera- 
tures. 

The positions and structure of the peaks in dM/dH 
show the magnetic fields at which there are maxima in 
the magnetization relaxation rate at a given temperature, 
applied field and magnetization. The dashed lines mark 
the positions of the three observed peaks showing their 
shift with temperature. The peak that occurs at the 
lowest magnetic field shifts gradually from H — 2.41 T 
to H = 2.48 T (by 0.17£?o) as the temperature decreases 
from 1.55 K to 1.23 K, that can be identified as the TAT 
regime. 

As the temperature is reduced from 1.23 K to 1.15 K, 
this maximum is shifted abruptly to H = 2.62 T (i.e., 
by 0.347?o) and remains approximately at the same value 
upon further decrease of temperature to 0.62 K. The mid- 
dle peak also shifts gradually from 2.81 T to 2.94 T in 
the interval between 1.55 K to 1.10 K, but it then shifts 
abruptly between 1.10 K and 0.91 K, at lower tempera- 
ture than the previous maximum. Below 0.91 K the po- 
sition of this peak also remains constant. We attribute 
this temperature independent regime to the pure QT. As- 
suming that at the lowest temperature only the ground 
state m esc — 10 participates in tunneling, we can iden- 
tify the positions of the observed peaks by the resonance 
field values H(n,m csc — 10), n = 5, 6, and 7, calculated 
according to Eq. (2), from low to high magnetic field. 
The peaks labeled n = 5 and 6 show an abrupt crossover 
between TAT to QT, which occurs in an interval of ap- 
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FIG. 2. Field derivative of normalized magnetization 
d(M /M s ) I 'dH vs H z at different temperatures measured at 
9 — 10° between the applied field and the easy axis of the 
crystal. The curves are offset for clarity. The dashed lines 
mark the positions of the maxima in dM/dH. Arrows corre- 
spond to the resonance fields H(n,m csc — 10) for n = 5, 6 
and 7 calculated according to Eq. (2). 
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FIG. 3. Peak positions (in the units of Bo = 0.42 T) versus 
temperature for 9 = 0° (squares), 9 = 10° (triangles down), 
9 = 20° (triangles up), 9 = 35° (circles). The bars on the left 
hand side of the graph show the escape levels calculated using 
Eq. (2). The accuracy with which the peak positions can be 
determined is approximately the size of the symbol. 



proximately 0.1 K. In contrast to n = 5 and n = 6 peaks, 
the n = 7 peak shifts to higher field (from 3.20 T to 
3.52 T) step-wise in the whole studied temperature in- 
terval (~ 1 K). We consider the crossover to QT in this 
case as gradual (second-order). 

Peak positions as a function of temperature, ob- 
tained from hysteresis experiments performed as de- 
scribed above, are summarized in Fig. ^. This plot shows 
the values of the longitudinal field, at which the maxima 
in dM/dH occur, versus temperature for the four stud- 
ied angles, 6 = 0°, 10°, 20°, and 35°. The peak positions 
were corrected for the effects of the internal field as de- 
scribed in Ref. |2^] . The bars on the left hand side of the 
figure show the escape levels calculated using Eq. (2), 
with parameters from spectroscopic data • 

Analyzing this graph, we can make following observa- 
tions. First, for larger angles, and therefore higher trans- 
verse field, peaks with lower indices (i.e., lower H z ) can 
be observed in the experimental time window. The low- 
est step observed for 6 = 0° and 10° is n = 5, for 6 = 20° 
it is n = 4, for 6 = 35° it is n = 3. This is consistent 
with the idea that the transverse field promotes tunneling 
and lowers the effective anisotropy barrier. We find that 
there is greater amplitude in lower lying peaks as the 
transverse field is increased. Second, two regimes can 
be distinguished: the high temperature regime, where 



the peaks gradually shift to higher fields with decreasing 
temperature, and the low temperature regime, where the 
peak positions are constant. We associate the first regime 
with the TAT and the second with pure QT. Third, the 
form of the crossover between these two regimes depends 
on the longitudinal field. For each sample orientation, 
peaks with lower indices show a more abrupt crossover 
between TAT and QT than peaks with higher indices 
(compare peaks n — 6 and n — 7 for 6 = 0° and 10°, or 
n = 4 and 5 for = 20°, or n = 3 and 4 for 9 = 20°). 

The experimentally determined escape levels can be 
compared to models [ ^3|JT^ , |2^ ] . As mentioned above, 
the relaxation rate can be written as a product of 
the tunneling rate between the pair of levels i and j, 
split by the energy Ay, and the Boltzmann factor: 
r y ~ A| eM-Eij/T), where E %] = {\E t + Ej\)/2 - E Q 
is the average energy of levels i and j above the low- 
est lying state in the metastable well Eq. The energy 
level m csc which dominates tunneling is determined as 
the level for which the tunneling rate T has a maximum 
at a given temperature, longitudinal and transverse field. 
The tunnel splittings Ay- were obtained by numerical di- 
agonalization of spin Hamiltonian (1), whose eigenvalues 
were determined with high precision p3fl . In Eq. (1) the 
second order uniaxial anisotropy term (i.e., BS Z ) was ig- 
nored for simplicity and a higher value of D (D = 0.66 K) 
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FIG. 4. Experimental escape levels m eS c as a function 
of temperature for n = 4 at two different orientations: 
a) 6> = 20° (triangles), and b) 9 = 35° (circles). Solid 
lines show m ea c calculated by numerical diagonalization of 
the Hamiltonian (1) with D = 0.66 K, B = K. 



FIG. 5. Minor loop hysteresis experiments for 9 — 10°, 
n — 6: a) typical hysteresis curve at T = 1.01 K; b), c), d) 
normalized d{M/M s )/dH at T = 1.01 K, 1.12 K, and 1.23 K 
respectively. The direction and order of the field sweeps is 
indicated by numbered arrows. 



was used so that the low lying energy level spacing and 
overall barrier height is better approximated. Figure 
|] shows experimental positions of the n = 4 peak at 
two different angles, 9 = 20° (Fig. 4(a)) and 9 = 35° 
(Fig. 4(b)) compared to the calculated m csc . Calcula- 
tions were performed using a Mathematica program for 
constant values of the longitudinal and transverse mag- 
netic field corresponding to n = 4, 9 = 20° and 9 = 35°. 

First consider the 9 = 20° data. A few major discrep- 
ancies between the calculated and experimental results 
can be noticed. The experimental crossover temperature 
is higher than the calculated one (T CX p = 1.25 K versus 
Tth = 0.7 K). Also the experimental temperature inter- 
val in which the crossover occurs is larger than predicted 
by the model. Experimental escape levels involved in the 
crossover (m esc = 7 — 10) are lying lower in the poten- 
tial well than the calculated ones (m csc = 5 — 10). The 
observed escape levels change more gradually than pre- 
dicted by the model: our observations show that none 
of the levels m csc = 7 — 9 are skipped in the crossover, 
while the calculations show that as many as three lev- 
els (m csc = 7,8,9) do not contribute significantly to 
the relaxation. Based on this comparison we conclude 
that there are additional mechanisms, such as transverse 



anisotropy, that enhance the tunneling rate of low lying 
levels. For 9 = 35° the crossover is more gradual than 
for 9 = 20° , in qualitative agreement with the theoretical 
model. As in the 9 = 20° case, experimental data show 
that lower lying levels (8 and 9) are active over larger in- 
tervals of temperature than predicted by the model. The 
overall better agreement with the calculated data in this 
case suggests that a large enough applied transverse field 
can become a dominant factor in determining the tunnel 
splittings. 

It is important to note that in our hysteresis experi- 
ments, in which the field is increasing, relaxation will first 
occur from thermally excited states because these come 
into resonance first (Eq. (2)). Therefore our data on 
Figures |^ and |] somewhat underestimate the crossover 
temperature. 

To show this we have performed a series of minor loop 
hysteresis experiments near the crossover temperature. 
In these experiments the field was swept at a constant 
rate back and forth across an interval, where a peak in 
dM/dH was observed. If the field is ramped towards 
positive saturation, relaxation from the excited levels 
is emphasized, if the direction of the sweep is towards 
lower magnetization, the ground state tunneling is fa- 
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vored. The results of these experiments in which the 
field was swept across the n=6 peak five times at 6 = 10° 
are shown on Fig. [|. Figure 5(a) shows a typical hys- 
teresis curve at T — 1.01 K. Figures 5(b-d) show plots of 
dM/dH at three different temperatures. For T = 1.01 K 
(Fig. 5(b)) the maximum in dM/dH for all sweeps occurs 
at 3.06 T, which corresponds to QT from m csc = 10. At 
a higher temperature, T = 1.12 K (Fig. 5(c)) the max- 
imum of the first sweep is at a lower field, H — 2.93 T, 
which can be identified as a thermal channel (m osc = 9). 
For the sweeps 2-5 the maximum occurs at the QT 
position (3.06 T), which suggests that at this tempera- 
ture QT still dominates. At higher temperature (1.23 
K, Fig. 5(d)) almost all relaxation occurs via a thermal 
channel at H = 2.87 T (m osc = 8), which shows that 
the crossover to TAT has already taken place. These re- 
sults show that regular hysteresis experiments emphasize 
the tunneling from the excited states, and therefore the 
crossover temperature determined from such experiments 
is underestimated. 

In conclusion, we presented low temperature magnetic 
hysteresis studies of thermally assisted and pure quan- 
tum tunneling in Mni2. The crossover between these two 
regimes can be either abrupt or gradual. These types of 
the crossovers are distinguished by the temperature inter- 
val within which the escape energy levels shift from ther- 
mally excited levels to the lowest state in the metastablc 
well. Through the hysteresis experiments we were able 
to identify these energy levels and study their behavior 
as a function of temperature. Our studies indicate that 
the higher longitudinal and transverse fields make the 
crossover more gradual. Comparison of our experimental 
data and theoretical model suggests that in Mni2 there 
must be additional mechanisms which promote tunnel- 
ing, such as a transverse anisotropy, which is responsi- 
ble for broadening the crossover for smaller transverse 
fields and increasing the tunneling rate from low lying 
levels. These experiments establish a lower bound on the 
crossover temperature (Fig. 3), as the standard hystere- 
sis experiments favor relaxation from the excited levels. 

We are grateful to Dmitry Garanin for fruitful discus- 
sions and help with numerical calculations. 
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